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Abstract 


In  this  paper,  the  constituent,  operation 
principle  and  regulation  process  of  the  bipro¬ 
pellant  variable  thrust  rocket  engine  for  a 
variable  area  injector  that  uses  a  solenoid 
valve  control  are  introduced. 

According  to  the  constituent  and  practic¬ 
al  operation  process  of  such  an  engine  control 
system,  the  mathematical  models  for  the  main 
links  are  given,  by  transfer  function  method, 
the  dynamic  and  steady-state  performance  is  ana¬ 
lysed,  according  to  the  calculation  results, 
the  effect  on  dynamic  performance  and  steady 
error  for  the  engine  control  system  is  further 
described. 


ii 


I. 

THE  CONSTITUENT  AND  OPERATION  PRINCIPLE  OF  THE  ENGINE  CONTROL 
SYSTEM 

In  this  text  we  want  to  introduce  the  use  of  a  pilot  circuit 
that  possesses  a  specific  frequency  and  can  vary  pulsewidth,  con¬ 
trolling  two  high  performance  solenoid  valves,  passing  through 
liquid  flow  rate  variations  to  change  the  pressure  in  the  hydraul¬ 
ic  pressure  cavity,  making  the  injector  needle  valve  move  up  and 
down  to  change  the  injector  circulation  cross  section,  control  the 
propellant  flow  rate  variations,  and  bring  about  thrust  variations. 
The  control  system  is  primarily  composed  of  the  pilot  circuit,  the 
solenoid  valve,  and  the  thrust  chamber  that  possesses  the  variable 
cross  section  injector.  Their  interrelation  is  shown  in  Fig.  1. 

The  pilot  circuit  is  composed  of  the  operational  amplifier, 
the  pulse  oscillator,  the  resonant  delay  oscillator,  and  a  solenoid 
valve  drive  circuit.  The  operational  amplifier  conducts  alternate 
operations  through  the  control  voltage  signal  and  the  combustion 
chamber  pressure  feedback  signal,  and  sends  out  a  voltage  deviation 
signal.  The  pulse  oscillator  sends  out  a  voltage  pulse  signal  of 
a  specific  frequency  and  width,  and  goes  to  control  the  resonant 
delay  oscillator.  The  resonant  delay  oscillator  conducts  the  vari¬ 
able  width  pulse  operations  that  restrain  the  frequency  of  the  pulse 
oscillator  and  the  operational  amplifier  to  export  the  degree  of 


Fig.  1.  DIAGRAM  OF  THE  OPERATION  PRINCIPLE  OF  THE  ENGINE 
Key:  (a)  Control  voltage;  (b)  Resonant  delay  oscillator; 

(cj  Solenoid  valve  drive  circuit;  (d)  Filling  mechanism;  (e)  Oper¬ 
ational  amplifier;  (f)  Pulse  Vibrator;  (g)  Electro-ignitor;  (h) 
Resonant  delay  oscillator;  (i)‘  Solenoid  valve  drive  circuit;  (j) 
Solenoid  valve  intake;  (k)  Rubber  gasket;  (l)  Spring;  (m)  Solenoid 
valve  escape;  (n)  Pilot  circuit  section;  (o)  Pressure  capsule;  (p) 
Combustion  chamber;  (q)  Valve  needle. 


voltage  deviation  to  control  it.  This  variable  pulse  signal  of  a 
certain  frequency  goes  through  the  solenoid  valve  drive  circuit  to 
carry  out  amplification,  and  spur  on  solenoid  valve  pulse  operations. 


The  two  solenoid  valves  and  the  injector  hydraulic  pressure 
cavity  are  linked  together.  When  the  deviation  signal  of  the  pilot 
circuit  makes  a  positive  voltage  signal,  it  controls  the  solenoid 
valve  intake  operations.  Pressure  rises  in  the  hydraulic  pressure 
chamber,  and  pushes  the  injector  needle  valve  to  rise.  The  in¬ 
jector  cross  section  expands,  the  propellant  flow  rate  increases, 
and  the  thrust  rises.  Conversely,  when  the  deviation  signal  of  the 
pilot  circuit  makes  a  negative  voltage  signal,  it  controls  the  sol¬ 
enoid  valve  escape  operations,  ejects  the  liquid  in  the  hydraulic 
pressure  cavity,  and  causes  the  pressure  in  the  hydraulic  pressure 
to  drop.  The  injector  needle  valve  is  under  the  effect  of  the 
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spring  force,  the  injection  flow  area  reduces,  the  propellant  flow 
rate  reduces,  and  the  thrust  reduces. 


The  control  operation  of  the  injector  needle  valve  uses  a 
combustible,  goes  through  the  intake  solenoid  valve  to  enter  the 
hydraulic  pressure  chamber.  Because  the  controlling  liquid  flow  rate 
is  very  small,  it  goes  through  the  exhaust  solenoid  valve  and  direct¬ 
ly  ejects  into  the  atmosphere.  A  system  can  also  be  designed  that 
makes  the  controlling  fluid  return  to  the  combustible. 


Fig.  2.  STRUCTURAL  DIAGRAM  OF  THE  SOLENOID  VALVE 

Key:  (a)  Escape  orifice;  (b)  Needle  valve  seat;  (c)  Needle  valve; 

(d;  Spring;  (e)  Coil;  (f)  Intake  orifice. 


Thus  it  can  be  seen  that  this  kind  of  variable  thrust  rocket 
engine  system  changes  the  size  of  the  injector  flow  cross  section 
to  change  the  propellant  flow  rate,  and  obtain  the  objective  of 
thrust  control.  The  engine  propellant  supply  system  does  not  con¬ 
duct  regulation.  The  characteristics  of  its  program  are  that  the  con¬ 
trol  system  is  simple,  the  operation  is  dependable,  but  the  regulation 
precision  isn't  high  enough. 


The  pilot  circuit's  operational  parameter  is  as  follows: 

voltage  of  the  pilot  circuit's  power  supply  V  —28  CV] 
controlling  voltage  of  the  circuit  V, -0~5  CVj 
pulse  oscillator  operational  frequency  /  — 50~100  [HzJ 
Escape  voltage  of  the  feedback  pressure  capsule  Vf  —0^—5  [V] 
Escape  voltage  of  the  operational  amplifier 

VA  -0.1-1. 1-2.1  CV] 
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(g) 


Fig.  3.  DIAGRAM  OF  THE  CONTROL  SYSTEM 

Key:(a)  Intake  solenoid  valve;  (b)  Control  voltage;  (c)  Control 
circuit;  (d)  Injector;  (e)  Combustion  chamber;  (f)  Exhaust  sol¬ 
enoid  valve;  (g;  Pressure  capsule. 

II. 

MATHEMATIC  EQUATIONS  FOR  THE  MAIN  LINKS  OF  THE  CONTROL  SYSTEM 

The  solenoid  valve  of  the  control  system  is  controlled  by 
the  pilot  circuit,  and  conducts  variable  width  pulse  operations. 
Therefore,  a  mathematical  model  of  the  pilot  circuit  and  solenoid 
valve  can  be  established.  In  the  mathematical  treatment  process 
we  must  encounter  the  nonconsecutive  alternate  system,  and  approach 
the  problem  of  the  Z  alternate. 

But  the  input  signal  of  the  control  system-control  voltage  Vc, 
and  the  output  signal- thrust  chamber  pressure  Pc  are  consecutive 
signals,  only  the  middle  links  are  pulse  signals.  In  order  to  sim¬ 
plify  the  overelaborate  mathematical  deductions,  we  can  build  the 
mathematical  models  of  the  main  links  of  the  circuit,  solenoid  valve, 
injector,  combustion  chamber,  etc.  based  on  the  actual  operations 
of  the  engine  system  and  the  regulation  conditions,  in  accordance 
with  consecuti /s  system  mathematical  treatments. 

1 .  Equation  for  the  Pilot  Circuit 

Basic  Assumptions: 

1)  The  pulse  operation  frequency  of  the  pilot  circuit  is  specific. 

The  size  of  the  circuit  operational  frequency  influences  the  dynamic 
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state  of  the  engine  system  and  the  steady-state  performance.  More¬ 
over,  the  size  of  the  frequency  receives  influence  from  the  per¬ 
formance  of  the  solenoid  valve  itself. 

Regarding  the  explanation  of  the  engine  system,  experiments 
show  clearly  that  its  operational  frequency  can  choose  within  the 
scope  of  50-100  Hz. 

2)  The  pulsewidth  of  the  pilot  circuit  is  in  a  direct  ratio  with 
the  voltage  deviation  (control  voltage  and  feedback  voltage  dif- 
ferences-deviation  signal).  (Experiments  prove  the  variations  of 
pulsewidth  in  direct  ratio  with  basic  voltage  deviation.) 

In  order  to  enhance  the  dynamic  state  performance  and  the  steady 
state  performance  of  the  engine  system,  the  pulsewidth  should  adjust 
as  wide  as  possible.  The  pulse  frequency  and  pulsewidth  of  the  pi¬ 
lot  circuit  do  not  affect  each  other,  and  can  be  separated  from 
regulation. 

3)  The  duration  of  the  starting  foreward  edge  of  the  solenoid  valve 
is  specific.  Because  the  solenoid  valve  is  an  inductive  element, 
the  solenoid  force  must  overcome  the  friction,  the  spring  force, 
the  damper  force,  and  the  hydraulic  pressure  of  the  needle  valve. 

The  solenoid  force  is  in  direct  ratio  with  the  amount  of  current. 

The  current  rises  to  the  duration  variations  in  accordance  with  the 
index  curve.  Therefore  it  forms  the  foreward  edge  when  the  solenoid 
valve  is  open.  The  curve  of  a  solenoid  valve  operational  current 
shown  in  Fig.  4  expresses 

The  starting  foreward  edge  =3-5ms 

The  needle  valve  opening  period  t/ =0 . 5ms 
n  represents  one  effective  pulsewidth  of  the  solenoid  valve  pulse. 

To  increase  the  pulsewidth  actually  is  to  increase  the  effective 
pulsewidth. 


Fig.  4.  CURVE  OF  A  SOLENOID  VALVE  OPERATIONAL  CURRENT 
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On  the  basis  of  the  above  analysis,  we  can  express  the 
equation  for  the  circuit-solenoid  valve: 


k- 

Vr 

o 


Where  is  circuit  gain 


(D 


[S/v] 


Eri=Vf>-V^  Differences  of  control  voltage  and  feedback  voltage  (Vj 
*  o 

n  is  effective  pulsewidth  jS^ 


When  V.=0,  E  =V 
X  o  c 


When  engine  regulation  arrives  at  operational  conditions,  after 
pulse  frequency  is  chosen,  the  pulsewidth  regulation  is  the  widest, 
also  then  is  set.  For  example,  f=50Hz,  Eq=Vc=5V,  to  the  pulse¬ 
width  that  should  be  the  largest  ,*•+ T*  =20ms ,  (  *•  =5ms ,  n=1  5ms ,  omit 
*'),  then  Kfe=  n/EQ=3x10"3  [S/Vj. 

The  equation  for  the  pulsewidth  of  the  solenoid  valve  is: 

b(t)«kb.E(t)  (2) 

The  effective  operational  duration  that  is  within  the  duration  of 
the  solenoid  valve  unit  is  called  the  effective  pulse  operational 
coefficient  of  the  solenoid  valve,  using  ?j(t)  to  express  clearly: 

u(t)=fb(t)  (3) 

Where  f  is  the  pulse  operational  frequency  of  the  solenoid  valve (HzJ 
V*(t)  is  the  effective  pulse  operation  coefficientC.’Kt )<1 ) 


2.  Equation  for  the  Flow  Rate  of  the  Solenoid  Valve 


Basic  Assumptions: 

1 )  After  the  solenoid  valve  goes  through  the  foreward  edge  interval 

within  the  pulse  operation  period,  the  needle  valve  opens  instan¬ 
taneously.  Because  after  the  solenoid  valve  goes  through  the  fore¬ 
ward  edge  interval  *vs3~5ms,  the  duration  that  the  needle  valve  is 
open  is  only  «r=0. 5ms, */«**#  hence  can  be  dismissed. 

2)  Experiments  show  that  the  closing  time  of  the  solenoid  valve 
needle  valve- the  back  edge  lessens  to  0.5ms,  can  be  omitted,  and 
deems  the  nee  le  valve  to  close  instantaneously. 


On  the  basis  of  the  preceding  assumptions,  considering  the 
structure  and  dimensions  of  the  solenoid  valve.  We  can  list  the 
equation  for  the  flow  rate  of  the  intake  solenoid  valve  and  the 
exhaust  solenoid  valve. 

The  equation  for  the  flow  rate  of  the  intake  solenoid  valve: 

fyzp.-pMWi 

Where  n  is  the  flow  rate  coefficient 

A  is  the  circulation  cross  section  of  the  solenoid  valve 
By  the  calculation  of  geometric  dimensions  of  the  solenoid  valve 
(reference  material ) , 

A-Z.llXD-i/Z') 

■* 

is  the  largest  interval  of  the  solenoid  valve  opening 
d  =0. 1/^0. 4mm 

D  is  diameter  of  thesolenoid  valve  escape  orificejjnm J 
Pq  is  hydraulic  pressure  of  the  solenoid  valve  intake 
orifice  £jkg/mm 

^H(t)  Pressure  °f*  the  main  hydraulic  pressure  cavity  of 
the  thrust  chamber  Qkg/mm^J 

The  equation  for  the  flow  rate  of  the  exhaust  solenoid  valve: 

The  liquid  of  the  exhaust  solenoid  valve  will  directly  enter  the 
atmosphere . 

The  solenoid  valve  is  pulse  operational.  Therefore,  the 
actual  flow  rate  of  the  solenoid  valve  should  lead  into  the  ef¬ 
fective  pulse  operation  coefficient  n(t),  and  then  the  equation  for 
the  flow  rate  of  the  intake  solenoid  valve  and  exhaust  solenoid 
valve  is 
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Ql0)-r,0)nAiJyLzPt-PH<.t')l 

Suostitube  the  above  equations  with-n>(t)  and  b(t)  of  equations 
(2)  and  (3)*  and  it  makes: 

(4) 

Then  the  preceding  equation  for  flow  rate  becomes 

o.co  *  KqEVWP'-PhO) 

Qt(t)-K9E0W'Pni0 

With  the  preceding  principle  knowledge  of  the  engine  system 
operations,  when  there  is  deviation  E(t)>0  of  the  control  voltage 
and  feedback  voltage,  it  controls  the  intake  solenoid  valve  oper¬ 
ations.  Conversely,  E(t)<0  operates  .the  exhaust  solenoid  valve. 
Which  one  operates  depends  on  whether  the  deviation  signal  E(t)  is 
positive  or  negative. 

In  order  to  use  an  equation  for  the  pilot  return  circuit  to 
express  the  equation  for  the  flow  rate  of  the  solenoid  valve,  we 
must  cite  two  functions.  Making 


Fig.  5.  GEOMETRIC  DIMENSIONS  OF  THE  SOLENOID  VALVE 
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‘<*>-2T£WC|£(')I+£(')] 


_ 1_ 

w(«)=2iiicoi  :i£(o i-^(03 


(6) 


Clearly, 


£(/)>0,  u(*)**l,  w(«)“0i 

E(t)<  0.  u(e)=*0,  — 


Then  the  equation  for  the  flow  rate  of  the  solenoid  valve  becomes 

<?,(0  •KQu(e)EOWP*-PH(t)  (7) 

Q,(0-tfow(«)£(OVTOI  (8) 


These  are  two  primary  non-linear  equations.  On  the  basis  of 
applied  design  requirements  of  the  engine,  each  time  the  extent  of 
thrust  regulation  is  comparatively  small,  thrust  regulation  vari¬ 
ations  are  a  smooth  transition.  We  can  use  extreme  progression 
to  carry  out  linearization  with  two  equations.  As  to  the  intake 
solenoid  valve: 


Where 


in  which 


■  gj$|( 

&*&&=*?•*** 


(9) 


.  fc. 


Mo) 
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makes 


(11) 


thus 


,,  jKZrKoK: 

Ki^Kikt 

QM-K'.EiO-ktAPM 


(12) 

(13) 


The  same  reasoning  obtains  the  equation  for  the  flow,  rate  of  the 
exhaust  solenoid  valve: 


where 


0  -  K  #(#) + K  tAPiit) 

jr;^nr«jr; 

'  K'iMfU'FZ 
Kt  -i.  ]-?■ 


(U) 

(15) 

(16) 

(17) 

(18) 


Equations  (13)  and  (14)  can  be  expressed  as  one  equation: 

Q(0-ol(0-o.(0-c^:-A,:)E(o-(/Cf+^:)APi,(o 


making 


k,~k:-k: 


KwKl+Kl 


thus 


Q(t)  -  K  .ECO 


(19) 


Clearly, 


when 

exo>o  K  *  ,0(f)-6iC0. 


EW<O.K.—Kl.K,m-Ki, 

Q(#)--Q.(  0 


The  corresponding  alternate  equation  for  the  flow  rate  of  the  sol¬ 
enoid  valve  is 


(20) 


10 


3.  Equation  for  the  Main  Capsule  of  the  Thrust  Chamber 


The  pressure  of  the  main  capsule  of  the  thrust  chamber,  ac- 
cording  to  design  requirements  P^( t )<20kg/cm  ,  and  the  comparison 
of  common  hydraulic  pressure  systems,  the  capsule  pressure  is 
comparatively  low.  In  the  equation  for  the  capsule  we  don't  need 
to  consider  the  distortion  of  the  capsule  housing  and  the  capsule 
climate  due  to  pressure  rising. 

Because  the  spring  has  preset  tension  Fq,  when  Fo  AHPHo' 
the  injector  needle. valve  starts  to  move.  Therefore  when  PH(t)=P 
y(t)=0.  Therefore 

Agdy(t)=Q(t)dt 

When  dy(t)>0,Q(t)=Q1  (t)  ;  dy(t K0,Q(t )=Q2(t) . 

y(t>=]4JQ(t)dt  (21) 

The  corresponding  alternate  formula: 

Y ( s )  =ip  . «tt Q ( s )  (22) 

AH  b 


Fig.  6. 


HYDRAULIC  PRESSURE  CAPSULE 
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4.  Equation  for  the  Injector  Needle  Valve  Movement 


A«P«t*> 


Fig.  7.  DIAGRAM  OF 
INJECTOR  NEEDLE  VALVE 
TENSION 


The  conditions  of  the  injector  needle  valve  tension  as 
shown  in  Fig.  7  are  expressed: 

FQ-Spring  preset  tension(kg) 

K-Spring  rigidity  [kg/mm) 
y(t)-Needle  valve  displacement Jmm] 

dy(t)/dt-Needle  valve  certain  moment  of  movement  speed  [mm/3 

2 

PjjQ-Hydraulic  pressure  when  the  needle  valve  starts  to  rise  (kg/mm  j 

APir(t) -Hydraulic  pressure  increments  £kg/mm^7 

“  "  2> 

Ag-Hydraulic  pressure  function  area  [mm  J 

B-Needle  valve  movement  damper  coefficient  Jkg»  s/mm"] 

Ac-Operational  area  of  combustion  to  the  needle  valvefinm  J 

2  ^ 
Pc(t)-Combustion  pressure  [kg/mm  J 

m-Needle  valve  mass jkg* s^/mmf 


The  equation  for  the  injector  needle  valve  becomes  (reference 
material^lj  ) : 


(23) 


The  corresponding  alternate  formula  is: 

(mS'+BS  +  JOY(s)-A.P.(s)+4HPM(») 


(24) 
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5.  Equation  for  the  Flow  Rate  of  the  Injector  Needle  Valve 


G.(0-^.^25vAPT  oxidant 
GfiO-ii/^rJZgYAPt  combustible 

Through  many  experiments  it  has  been  proved  that  when  the 
displacement  of  the  injector  needle  valve  is  comparatively  small, 
the  changes  of  flow  rate  and  displacement  cause  a  linear  relation¬ 
ship.  The  two  flow  rate  equations  can  be  expressed  as 


<?.(  0-#.y(0 

(25) 

GK  0-#,y(0 

(26) 

In  the  formula  the  coefficients#* and #, can  traverse  the  experimental 
determinant  of  the  injector  flow  rate. 

The  corresponding  alternate  formulas  for  the  flow  rate  equations: 


G.<S>-#*K(S) 

(27) 

G/(S)-#/KCS) 

(28) 

6.  Equation  for  the  Combustion  Chamber 
Basic  Assumptions: 

1)  Combustion  delayr  is  constant.  Consider  that  after  the  propellant 
is  injected  into  the  combustion  chamber,  it  goes  through  the r period, 
and  Combustible  gas  is  instantaneously  produced. 

2)  Combustible  gas  is  the  ideal  gas. 

3)  The  speed  of  the  flow  rate  variation  and  the  speed  of  pressure 
wave  propagation  are  comparable.  Consider  that  the  speed  of  pressure 
wave  propagation  is  infinitely  great,  and  the  pressure  in  the  com- 


bustion  chamber  in  any  moment  or  place  is  even 


On  the  basis  of  the  above-mentioned  assumptions,  we  can  ob¬ 
tain  the  equation  for  the  combustion  chamber  (reference  material  1  ) 


jjfr  *0*  ° + XP<(0  -  ~  + G/(*  -t0  , 


(29) 


Where  -Combustion  chamber  volume  [inm^J 
?.-Ideal  gas  constant  (mm/°K| 
r#-Combustible  fuel  temperature  C*K] 

2 

^4,*-Injector  nozzle  cross  section  area  [mm  J 
0 -Combustion  chamber  comprehensive  parameter's] 


The  corresponding  alternate  formula  : 


7.  Equation  for  Feedback  Pressure  Capsule  and  Circuit  Operation 
Amplifier 

The  pressure  capsule  of  the  engine  system  being  used  is  a 
resistance  pressure  capsule.  For  this  kind  of  pressure  capsule  we 
may  consider  it  an  amplify  link.  Its  equation  is: 


K/(f)  -K(P.OY 


The  equation  for  the  circuit  operation  amplifier  is: 


(31 ) 


F(0-K.(0-K/(0 


The  corresponding  alternate  formulas  are: 


Vii^-KtPM 


.tfCO-P.CO-KKO 


(33) 


(34) 
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TRANSFER  FUNCTION  OF  THE  ENGINE  CONTROL  SYSTEM 

To  synthesize  the  preceding  equations  (20),  (22),  (24),  (27) 
(28),  (30),  (33),  and  (34).  the  diagram  shows  the  interrelation 
between  several  variables. 


Fig.  8.  DIAGRAM  OF  THE  CONTROL  SYSTEM 


In  which  j=vc/RTc,  q~A,/0 

Because  the  engine  uses  bipropellant  spontaneous  combustion 
thrust,  ignition  delay  timer  is  comparatively  small.  General 
?<5ms  can  be  considered  e**  =  1  . 

The  volume  of  a  small  engine  combustion  chamber  v,  is  compara¬ 
tively  small,  and  the  temper.ature  t,  is  comparatively  high.  RT,>v, 

can  deem  j  =  «./RT  =0. 

c 


Fig.  9.  SIMPLIFIED  DIAGRAM 


The  simplified  diagram  is  the  following  formula: 


where 

K.— k#.+  #/)--f<#.+#/) 

q  A, 

Kc  represents  combustion  chamber  gain  [kg/mm^] 

The  close  link  transfer  function  of  the  control  system: 

^(,)  “«j* + +k.K.K  ,AHIK~r» 


(35) 


(36) 


(37) 


2 

m-Injector  needle  valve  mass  kg  s  /mm 
B-Injector  needle  valve  damper  coefficient  kg  s/mm 

2 

Ag-Needle  valve  hydraulic  pressure  operation  area  mm 
K^M-Solenoid  valve  pressure  coefficient  mm  Vs  kg 
-Solenoid  valve  flow  rate  gain  mnr/s  V 
K-Spring  rigidity  kg/mm 
K,  -Combustion  chamber  gain  kg/mm 

2 

A,  -Needle  valve  fuel  operation  area  mm 

-  2 

-Feedback  pressure  chamber  gain  mm  V/kg 

When  calculating  the  transfer  function,  when  using  B,Pjjt(t=0), 
^Ht^Ho  and  Kc*  according  to  the  actual  structure  of  the  engine 
and  the  operational  parameters,  we  can  use  experimental  method 
determinants, 

IV. 


ANALYSIS  OF  THE  ENGINE  CONTROL  SYSTEM  PERFORMANCE 


1 .  Analysis  of  Steady  State  of  the  Control  System 

Characteristics  of  the  control  system  transfer  function: 
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+  ( B  +  A;./ .V ,J(  )S  +  K-K  .A.  +  A  .A ,K /A*/K  ftt  -  0 
I'-’B  +  Ait/K  ,w,v~A  -k.At  +  K'KtkfAii/K 


Root  of  the  characteristic  equation: 


S.- 

— A  +  v  A*— 4«i* 

2m 

(38) 

s,- 

— A— VA*  — 4m» 

1m 

(39) 

Because  A>0,  discussing  the  stability  of  the  control  system 
we  only  must  discuss  the  root  of  the  characteristic  equation*  and 
that  is  enough.  From  the  root  of  the  characteristic  equation 
we  can  see  that: 

1)  v<0,  thus  there  is  the  solid  root  of  S^>0.  The  transition  weight 
of  output  monotonically  increases  according  to  the  exponential  law, 
and  the  system  is  unstable.  However,  from  the  plan  of  the  engine 
control  system  we  can  see  that  if  the  unstable  condition  emerges, 

the  pressure  in  the  combustion  chamber  monotonically  increase  or 
monotonically  decrease.  Because  the  control  system  has  a  closed  link 
system  of  feedback,  the  pressure  in  the  combustion  chamber  exceeds 
the  specific  value  made  to  correspond  with  the  solenoid  valve  opera¬ 
tions.  This  kind  of  combustion  chamber  pressure  will  emerge  accord¬ 
ing  to  specific  frequency  and  amplitude  oscillation  model  operations. 
The  exhaust  and  intake  solenoid  valve  alternate  operations. 

2)  v>0,  the  solid  root  of  A'  -4mv>0,  S^O,  the  control  system  is 
stable . 

3)  v>0,  A*-4mv<0,  S1  and  Sgboth  are  the  conjugate  compound  radical 
of  the  characteristic  equation. 

In  order  to  further  explain  that  these  parameters  are  the  pri¬ 
mary  influential  factors  of  engine  control  system  performance,  we 
will  use  the  parameters  of  the  engine  to  conduct  actual  calculations 
to  add  to  the  explanation. 
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Primary  parameters 


Table  1  . 


e 


Symbol 


Numerical  Value 


1  .2x10 


2 . 5x1  0 


1  .2x10 


When  experimenting  with  water  as  the  hydraulic  pressure  opera 
tive,  we  can  calculate  the  condition  of  the  thrust  rise  according 
to  the  data  given  in  Table  1 : 


t6.78x 

K\  -s/pt-P*,  -0.S2Ck*/8l«lI,3,', 
Kf-y  £'.CN/p,-Pw.  3",-0.a2CV-Knj/kc,M] 
.  K.-K0K:  -2.17  x  10*E*b//i  VJ 
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If  the  injector  design  is  good,  when  the  needle  valve  dis¬ 
places  h_ov<b.25mm,  the  variations  of  injector  flow  rate  and  needle 

ulaX 

valve  displacement  are  similar  to  a  linear  relationship.  The  com¬ 
bustion  chamber  gain  K±  can  be  expressed: 


If  Pc^—lOx  icr*kg/mat,  h  =0.25mm  thus  O.iOkg/mat*. 

max 

In  accordance  with  the  preceding  parameters  we  are  able  to 
calculate  the  values  of  the  roots  v  and  A  of  the  characteristic 
equation,  and  determine  the  stability  of  the  engine  control  system. 
Moreover,  we  can  see  the  influence  of  various  parameters  on  trans¬ 
fer  function,  and  further  analyse  the  primary  influential  factors  of 
the  performance  of  the  engine. 

*~K-K*d.+K.KJC,AM/K,H<ml20~M+Z5S.SSx  10*-35.578x  lO'kg/mm 
A- B +A'h/K,h  -3.8+ 11 . 797  X  10*  - 1 1 .83  x  10*kg •  s/min 
A*-4nw- 139.9489  X  104- 1 .71- 139.487  X  10*(kg-s/nim)* 


Calculating  the  expression  v>0,  A*-4mv>0,  the  root  character¬ 
istic  equation  of  the  transfer  function  is  a  negative  root.  Therefore 
the  engine  control  system  is  stable. 

2.  Analysis  of  the  Transision  Period  of  the  Engine  System 

Because  the  mass  of  the  engine  injector  needle  valve  is  Vo 
small,  m4<B=Ag/Kp£,  we  can  omit  the  quadradic  sum.  The  pressure 
variations  of  the  combustion  chamber  that  the  temperature  testing 
instrument  measures  monotonically  rise  and  fall  according  to  the 
exponential  curve.  Because  of  this,  the  engine  control  system  can  be 
simplified,  and  the  transfer  function  formula  changes  to: 


40) 


KJCMKpu-JL 

(A/»)S+i 


9>*KtK'AH/KpM’*  is  the  amplitude  coefficient 
r-A/*  is  the  duration  constant 

According  to  the  analysis  of  a  link,  when  the  combustion  chamber  pres¬ 
sure  p,  arrives  at  the  95 %  of  specified  value,  its  transistion 
period  tp—ZT,  Therefore, 

K - (40) 


Because  the  variable  scope  of  one  parameter  in  itself  is  not  large, 
from  the  calculations  of  numeric  value  of  the  front  side  we  see  that: 

\K-K*A.\4:K.KJC/Am/K,m 

The  interrelated  expression  of  the  transition  period  *p  can  be 
simplified  as: 


tf-K3ZJT/  (41) 

By  equation  (41)  we  can  see  that  the  transitional  period  is  primarily 
dependent  on  the  injector  needle  valve  hydraulic  pressure  operation 
area  Ag  and  the  solenoid  valve  flow  rate  gain  In  order  to  re¬ 

duce  the  transitional  period  ip, ,  the  injector  needle  valve  hydraulic 
pressure  operation  area  Ag  must  fully  decrease.  The  properly  en¬ 
larged  solenoid  valve  circulation  cross  section  A,  increases  the 
pulsewidth  or  raises  the  pulse  operational  frequency.  Both  are 
effective  methods  for  reducing  the  transition  period. 

3.  Analysis  of  the  Engine  System  Steady  State  Error 

By  Pig.  9,  the  diagram  of  the  engine  control  system,  we  can 
arrive  at  the  formula  for  the  deviation  signal  E(t): 
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where 


E(*> — 


K.K-AulK .. 

*<•>  -  **.+  is +£/£">  +z=zx 


In  the  engine  control  system  that  this  text  explains,  control 
voltage  V,(t)  intakes  by  the  form  of  rising  function.  Therefore, 
here  we  will  only  discuss  the  seat  deviation  of  the  rising  intake 
signal.  V/e  will  not  grant  discussion  to  the  velocity  error  and 
acceleration  error. 


where  K r  is  control  voltage  amplitude  coefficient. 

Using  the  alternate  whole  value  theorem  we  can  arrive  at  the  common 
expression  for  steady-state  error: 


jr 


in  - 


Substituting  W(o)  we  can  arrive  at  the  expression  for  the  engine 
control  system  steady-state  error: 


In  the  same  way,  jlT— K^i\4C.KiK',KtA*IKrm 

The  steady-state . error  expression  can  be  simplified  to  the  following 
form: 


M-Txj&rffzr,*’' 

From  equations  (9)»  (10),  (11),  and  (12)  we  obtain: 
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Thus  the  expression  for  steady-state  error  is: 


KO— 2  KJi  /AhIp'-^ 


Increase  the  spring  rigidity  K,  and  reduce  the  injector  needle 
valve  hydraulic  pressure  operation  area,  and  both  can  cause  the 
steady-state  error  l(t)  to  increase. 

The  thrust  regulation  weight  is  great,  the  deviation  signal  Et 
is  large,  therefore  steady-state  error  is  great.  There  is  more 
steady-state  error  when  there  is  high  thrust  (P^  is  more)  regula¬ 
tion  than  when  there  is  low  thrust  (P^  is  less)  regulation. 

The  preceding  discussion  is  aimed  at  the  condition  of  thrust 
rise.  As  to  the  condition  of  thrust  decline,  solenoid  valve  flow 
rate  gain  Ke  and  solenoid  valve  pressure  coefficient  Kpp  use  related 
equations  (15),  (16),  (17)  and  (18).  If  we  use  the  parameters  of 
the  exhaust  solenoid  valve,  other  parameters  do  not  change,  and  we  can 
obtain  useful  analytical  results. 

We  still  must  point  out,  regarding  the  analysis  for  regulation 
performance  of  an  engine  system,  we  have  not  considered  the  steady- 
state  error  of  the  control  system  element  in  itself  due  to  variations 
of  structure  and  parameters  influencing  the  rising  of  its  degree  of 
sensitivity.  For  instance,  the  solenoid  valve  needle  valve  movement 
intervals  increase,  the  magnetic  block  increases,  the  foreward  edge 
of  the  needle  valve  opening  enlarges,  and  causes  steady-state  error 
to  increase.  In  the  same  way,  the  solenoid  valve  intake  pressure  Po 
rises,  the  starting  force  increases,  and  can  also  cause  the  starting 
forward  edge  of  the  solenoid  valve  to  enlarge,  influencing  the  per¬ 
formance  of  the  dynamic  state  and  the  steady-state.  In  addition,  the 
non-linearity  of  the  feedback  pressure  capsule  can  also  influence 
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the  regulation  performance  of  the  engine  control  system.  Similar 
factors  have  not  been  included  in  the  above  analysis. 
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